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ABSTRACT: 


The radiative lifetime t for the decay of massious neutrinos is 
calculated using various physical models for neutrino decay. The results are 
then related to the astrophysical problem of the detectability of the decay 
photons from cosmic neutrinos. Conversely, the astrophysical data are used to 
place lower limits on t. These limits are all well below predicted values. 
However, an observed feature at - 1700 ^ in the ultraviolet background 
radiation at high galactic latitudes may be from the decay of neutrinos with 
mass - 14 eV. This would require a decay rate much larger than the 
predictions of "standard" models but could be indicative of a decay rate 
possible 'n composite models or other new physics. We may thus have found an 
important test for substructure in leptons and quarks or other physics b«»yond 
the standard electroweak model. 

Subject headings: cosmology - elementary particles - neutrinos - 

ultraviolet: spectra 


PRECEDING PAGE BLANK NOT RLMED 


3 


I. INTRODUCTION 

While suggestions tying astrophysical observations with the possibility 
of massious^ neutrinos have been around for some time (Gerstein and Zel'dovich 


1. The literature heretofore on this subject uses either the phrase "non-zero 
neutrino rest mass" or the adjective "massive". The former phrase is 
technically accurate but cumbersome. The adjective "massive" in standard 
English usage connotes heavy, bulky, or large objects, hardly appropriate for 
m^ s adjective "massy" seldom used, has an even more pronounced 

connotation in this direction. We prefer here to coin the adjective 
"massious" (having the property of mass) as more appropriate for present 
physics usage. 


1966, Cowsik and McClellend 1972, Szalay and Marx 1976), the advent of grand 
unification theories (see, e.g., Langacker 1981) and (as we will suggest here) 
composite models of quarks and leptons (see review by Harari 1980) as well as 
recently reported experimental results implying finite (Reines, Sobel and 
Pasierb 1980) and cosmologically significant (l.yubimov, et al. 1980) neutrino 
masses, are stimulating much interest and work on the subject of massious 
neutrinos and their cosmological implications (Schranw and Steigman 1981, 
Dolgov and Zel'dovich 1981). We begin with a brief summary of the basic 
cosmological setting for a discussion of this topic. 


II. COSMOLOGICAL SETTING 

Since the radiative lifetimes of light massious neutrinos are expected to 
be much larger than the age of the universe, both from theoretical (de Rujula 
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and GlashoM 1980) and some observational (Cowsik 1980) considerations, one 
must look for the most copious source of neutrinos In the universe In order to 
look for photons from their decay. This source Is the big-bang Itself. For 
each neutrino flavor f and hellcity e^, the number density of neutrinos plus 
anti neutrinos In the universe Is 

n^ = 1.1 X 10 { 2 ^ ) cm (1) 

(see, e.g., Weinberg 1972). 

The total number density Is thus 

n = 110 zg, (2) 

V f 

taking T = 2.7K and the total mass is 

Denoting the fraction of the closure density of the universe 

In neutrinos, it follows that 

a = 0.01 h"^ z (4) 

V 0 V 

where h^ Is the present Hubble constant In units of 100 km s“^ Mpc"^ and 2^ Is 

In eV. Thus a value for 25 < z < 100 eV could close the universe 

V 

(0.5 ~ h^ ~ 1). We may compare equation (4) with the various values of tl 
associated with objects on different astronomical scales. The ratio of mass- 
to-llght based on dynamical mass measurements Increases with the Increasing 
scale size. It Is found that over distances much larger than typical 
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Interstellar scaUs, M/L Is proportional to scale size 
(M/L « r) up to distances of the order of ^ 1 Mpc (Davis, Tonry, Huchra and 
Latham 1980). Our version of Figure 2 of Davis et al . (1980), which takes 
account of additional data (Blitz 1979; Hoffman, Olsen and Sal peter 1980) Is 
shown In Figure 1. The curve shown In Figure 1 gives a functional 
approximation to the data of the form 

(M/L) » UQ[l-exp(-r/A)] (5) 

In solar units. At extragalactic distances, the h dependence Is also shown on 
the scale. The function (5) has the virtue that M/L « r for r << A and M/L ♦ 
const for r >> A as required by the observational constraint n ~ 2. The 
value for M/L corresponding to the critical density (I.e., fl » 1) Is shown by 
the circle marked C. It can be seen that there appears to be a scale size A ~ 
a few Mpc which Is characteristic of the non-lumlnous mass In the universe. 
This size Is Interestingly close to the galaxy clustering size ~ 4 Mpc 
(Peebles 1980) and Is of the order of the Jeans length (scaled to the present 
time) which one would obtain from the growth of gravitational perturbations of 
neutrinos In the range 

a few eV < m < a few tens of eV (6) 

~ V ~ 

(B1snovaty14:ogan and Novikov 1980; Bond, Efstathlou and Silk 1980; 
Doroshkevich, et al. 1980; Sata and Takahara 1981) This range of masses Is 
also relevant to the c^namlcal studies of Tremaine and Gunn (1979). It should 
also be noted that cosmological neutrinos can undergo violent relaxation 
(Lynden-Bell 1967; Steigman, Sarazin, ()u1ntana and Faulkner 1978) to produce a 
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m9 

density distribution n^ « as implied by rotation curve studies of the 
outer parts of galaxies (Blitz 1979) and that such a density distribution, 
when extrapolated to galaxy clusters, can give the observed relation 
M/L « r . It may also be noted that massious neutrinos ir the mass range (6), 
could close the universe (sec equation (4)) and thereby "solve” the "flatne's 
problem" as proposed by Guth (1981). (Of course, many workers do not share 
the view that fl must equal 1 and do not recognize a flatness problem). 

Without getting into such controversial areas as to whether or not a = 1 or 
whether neutrinos cluster on the scale of galaxy clusters, galaxy halos, or 
both, we will therefore concentrate our further discussion on the radiative 
decay of neutrinos in the mass range (6) and the consequences of searching for 
the decay photons. 

III. ASTROPHYSICAL NEUTRINO FLUXES AND RADIATIVE LIFETIMES 

It has been pointed out by De Rujula and Glashow (1980)that the 
wavelength range to search for photons from the decay of cosmological ly 
produced neutrinos (mass range given by (6)) lies in the far ultraviolet. 

This is because for the decay from a heavier (v') to a lighter mass (v) 
neutrino 


v' ♦ V + Y 


the emitted photon has an energy 



(7) 


( 8 ) 


in the rest system of the decaying neutrinos. 
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The neutrinos have been "adiabatlcally cooled" by the expansion of the 

universe so that their velocity spread is determined by the <1ynamics of their 

gravitational interaction rather than by thermal velocities. Typical 

velocities f ^ r neutrinos bound in galaxy halos would be ~ 300 km/s. For 

neutrinos in gala)^ clusters, the dynamical velocities would be ~ 10^ km/s. 

Thus, for Eq corresponding to a wavelength Xq ~ 1000 X (E^ ^ 12eV) the Doppler 

0 

spread of the lines would be ax ~ 1 A for neutrinos in galaxy halos and 
0 

~ 3 A for neutrinos in galaxy clusters (ax/x^ -x v /c). For the case where 
m' >> m, which might be expected in light of the large mass differences known 
to exist among the charged leptons, equation (8) reduces to 

E^ » m'/2, m' » m (9) 

In contrast to the narrow monochrotMtic radiation expected from nearby 
objects, there should also be continuum radiation at E < E^ (X > x„) from the 
decay of neutrinos which occurred in tne past when we integrate the line 
emission over all redshifts. 

The formulas for the as trophy si cal photon fluxes are as follows: 

1) The diffuse line intensity from the galactic halo is given by 

^X * t^ t ~ AX ~ ^ ^ 

where t and n' are the lifetime and density of v' neutrinos and the integral 
is along the line-of-sight of the telescope. 

2) The flux from an extragalactic source such as the halo of a 
nearby galaxy or a nearby galaxy cluster is given by 

F - — m / n'dV - — cm"^ s‘^ A 

4«R^tAX 4wR^tAX 


( 11 ) 
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where the volume Integral 
N = / n'dV 


( 12 ) 


gives the total number of v neutrinos In the source and R is the distance to 
the source. If the mass of a galaxy cluster or halo is assumed to be mainly 
from \>2 neutrinos, then 


2. X 

^ * STTeVT 


(13) 


where is the total mass of the source, usually given in solar mass (Mg) 
units. 

3) The continuum flux from the decay of cosmological neutrinos is 


1(E) 


c jh ^ «ni*z)c-EQ] 
’’**0 ’ 0 (1+zl (l+nz^'iz 


ci? s?! ejl 


(14) 


(see, e.g., Stecker, 1971) 

where z^. is the critical redshift of absorption of the UV flux. 

Since E » K/x where K ■ 1.24 x 10"^ eV X « he, in wavelength units and for 
Xq « hc/E^, equation (14) becomes 


I 


X 


cn Q 
4 iMqt 



[1 + (n-l)(l-XQ/x)] 



*0 ‘ “ \)'l * 'c* 


(15) 
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or, in numerical units (Stecker 1980; Kimble, Bowyer and Jacobsen 1981). 


3/2 


- 7.8 X 10^®hQ^ “ (n-l)(l-XQ/A)J“^^ 


cm s sr A 


( 16 ) 


Since the expected ultraviolet fluxes (10), (11) and (16) are 
proportional to the neutrino decay rate , the physics of neutrino 
decay (t for v* ♦ v + y) and the astrophysical observations are both related 
to the problem of determining the lifetime of putative massious neutri'ios in 
the mass range (6). 


IV. MODELS FOR RADIATIVE NEUTRINO DECAY 


To compute the radiative neutrino decay rate r ■ t~^, ive first note that 
the most general form for the amplitude is 


T(v‘ ♦ V ♦ y) • i ^ ? (P“Q) Oyy e%'*(a + bYg) T(p) (17) 


2 2 2 2 

where p - m' , (p-q) » m . a and b are dimensionless numbers while A 

characterizes the relevant mass scale, or combination of mass scales, involved 
in the decay interaction. Equation (17) follows fror' gauge and Lorentz 
invariance and leads to the decay rate 


r 


« (,^|2 ^ 
2A^ " 


(18) 


If m' » m. 
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r ■ 3.65 X 10“^^ (lai^ ♦ Ibi^) ev 


LA(6eV)J‘ 


(19) 


or 


T - 1.80 X 10® (lal^ + Ibl^)’^ sec. (20) 

Lm‘(eV)r 

Equation (20) Is the basis for a discussion about the lifetimes predicted In 
various models. The models have a wide range of characteristics, and It Is 
useful to characterize them by the parameters a, b, and A. 


A. Conservative GWS 

It may yet turn out that neutrinos really are massless and hence 
do not oscillate, as In the standard Glashow-Uelnberg-Salam (GWS) model with 
no right-handed neutrinos. In this case, a ■ b ■ 0, and 

T - - (GWS) (21) 

and there would be no more story to tell. This would also be the case for 
massfous neutrinos with conserved lepton number (flavor). 

B. Extended GWS 

On the other hand. It Is easy to extend GWS to Include neutrino 
masses and mixing. (The mass eigenstates differ from the weak-interactlon 
basis V|^.) Neutrino electromagnetic decay can now proceed by an Intermediate 
state consisting of a weak boson W and a charged lepton t both of which can 
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couple to the photon (see Figure 2). For three generations (1 ■ 1, 2, 3; t ■ 
e, Ml t) of Dirac neutrinos, and for m 3 > mj, for exan^le, 
the vj ♦ vj + Y decay rate Is (Marciano and Sanda 1977) 


r(v. 


Vj + y) 


9G^ a 
2048? 


2 2 

*’'3 " ***1 3 2 2 

(-~r — k ^ (m‘ ♦ mf) 

m 3 3 1 


2 2 2 2 2 2 2 
m^ s, ♦ ro c, - m^ m^ - m 

X [ , s^CjSj 


"T 

% 


( 22 ) 




In terms of the Kobayash1-Maskawa-1 Ike neutrino mixing 

angles (s^ = sin 8 ^i = cos 8 ^) (Kobayashi and Maskawa 1973) with no CP 

violation. For a general v* ♦ v + y, the scale Is 


A 


(Gpm* ) 


•1 _ 10 


14 


iirreVT 


GeV 


(23) 


and the numbers a and b are (Ignoring $^, c^ factors) 


lai, Ibl - - 4.3 X 10‘® (24) 

32i my 


This Is consistent with (22). Therefore we have 


T • — — i (Extended GWS) (25) 

[m‘(eV)r 

It must be remembered that the mixing angles may Increase this significantly 
(e.g., T ■ • If 8j ■ 0!). 


C. Heavy lepton 

The leptonic version of the Glashow-Il lopoulos-MalanI (GIM) supnresslon 
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mechanism was operative in (22) and led to the 0 numbers in (24). We 

can therefore achieve a larger decay rate by going to some model Involving 
heavier leptons. Fquatlon (24) Is then changed to 

lal. |b| - - 10"^ (26) 

32ir‘^ 


and so 


T 


CmMeV)]- 


(Heavy Lepton) 


(27) 


This agrees with detailed model calculations with an additional very heavy 
lepton (Pal and Wolfenstein 1981; Aliev and Vysotsky 1981) (fourth generation) 
and was first estimated by de Rujula and Glashow (1980). 


D. GIM-less 

Models where the GIM mechanism Is absent altogether could also decrease 
the lifetime to the order-of-magnitude (27) (de Rujula and Glashow 1980; Pal 
and Wolfenstein 1981/. We may write 

10^^ s 

T - --iii * (NO GIM) (28) 

[m'(eV)r 

with the caveat that this, as veil as our other estimates, could be signifi- 
cantly larger If mixing angles are sufficiently small. 

E. MaJorana«D1rac Neutrinos 

We may try to evade GIM suppression by considering _b^h Dirac and 
Majorana mass terms In the Lagranglan, a circumstance which can arise In 
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certain grand unified theories where the Majorana masses can be Induced by 
radiative corrections. Cheng and LI (1980) have studied the rates 
for u 6Y for these general neutrino mass eigenstates In an extended GWS 
model, and we can adopt their work to v* ♦ vy* If all six of the masses are 
small, we still have GIM cancellations. If we choose three of the masses to 
be as large as we wish, a fine tuning of the parameters In a most general mass 
matrix can enhance the decay rate. However, we still see the same lower limit 

37 

j I , — 5 (Majorana/Dirac) (29) 

[m' (eV)]^ 


F. Higgs 

Pal and Wolfenstein (1981) have also pointed out that Higgs Intermediate 

2 

states could enhance amplitudes by a factor of (M^/M^) where Is a Higgs 
mass. If there Is no GIM-llke cancellation In the remaining factors, then we 
can optimistically guess that 

T ~ (Higgs) (30) 

^ [m’(eV)]^ 

In the case where = 0.1, a four order-of-magnitude reduction would 

^ w 

result. 

G. Composite Models 

There has been much effort In recent years constructing composite models 
of quarks and leptons out of more basic particles. The research area Is quite 
new. We have no single calculation to offer as a good Indication of what to 
expect for a decay lifetime. However, we are able to obtain order of 
magnitude estimates for lifetimes If a scale A (composite size A*^) is 


given. This holds for a reasonably large class of models, an Important 
consideration since we want to be sure that there Is no general principle 
which states that the electromagnetic decay rate for composite neutrinos Is 
vanishingly small. 

The reasons for believing that the leptons and quarks are bound states o 
something else, have to do with the prol Iferatlon of particles and parameters 
with the fact that the higher generations resemble excitations of the "ground 
state" generation, with the mismatch between fundamental supersymmetry multl- 
plets, and known particles, and so forth (Hararl, 1980). Responding to such 
Incentive, many composite models have been proposed and, although neutrino 
radiative deca^y has not been studied, most of these models do embrace non-zer 
mass neutrinos and neutrino flavor nonconservation. The problems In bulldinc 
very light and very small particles with pointllke magnetic monents have not 
yet been solved, but only the scales for each class of models are needed for 
order-of-magnitude lifetime estimates. 

The composite scale ^ Is unknown, with the lower limit 

A > 1-10^ TeV (31) 

based on the absence of non-QED anomalous magnetic moments, on the absence of 
structure In scattering, and on Higgs compositeness. The limit from the 
absence of proton decay Is assumed to be no more severe (Chandra and Ray 
1980). The y ■» ey decay limit Implies A > 10® TeV In first-order and A > 10^ 
TeV In second-order (Barblerl et al . 1980). 

The standard weak-induced decay models have a suppression mechanism due 
to the conflict between the spin flip of the radiative decay and the chiral 
coupling of the weak boson. Therefore the transition magnetic moment fights 
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the V-A coupling (a coupling that Kwuld fix the helldty of a massless fermion 
line) so that the amplitude Is proportional to and vanishes with the neutrino 
masses. This Inhibitory factor Is separate from phase space, the magnetic 
dipole factor, and ar\y GIM suppression. The situation can be different In a 
composite theory. We still must have the overall spin flip, but the role of 
chiral symmetry Is less certain. If there Is no chiral symmetry Involved In 
the dynamics, this suppression mechanism Is not operative. Even If there are 
chiral symmetries of “naturalness", some of them must be broken by a set of 
nonvanishing mass scales. The constraint that fixes the composite wave 
function spin In a chiral limit may be satisfied by the presence of the other 
scale factors. 

If the dimensionless functions (of ar\y mass ratios In the dynamics), a 
and b In equation (18), entail chiral suppression, a and b are of the form 

a H ^ a’, b 3 ^ b' (32) 

with the mass scale A' possibly different from G. The radiative decay rate Is 
(from (18)) 

r - g, g = ^ (la|^ ♦ Ibl^). (33) 

A 

For the chiral case, we define 

G - (^)^ G', G' 3 (la'I^ ♦ Ib'l^) (34) 

We see that two mass powers come from the magnetic dipole q-f actor, two powers 
from the chiral suppression, and one power from phase space. The radiative 
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lifetime Is 


9.0 X 10^^ [A(TeV)]2 _ 

T ■ 7 > - " - ■■ ■ 5 

[m‘(eV)3^ 


( 35 ) 


or, from (34) 


9.0 X 10^^ [A(TeV)A'(TeV)]^ . 

^ ^ — [;;w 


(36) 


In the composite approach, protons are composites of composites and there 
are various ways In which Its decay may be Inhibited, with no direct Impli- 
cation for v' ♦ VY decays. On the other hand, u ♦ ey Is much more closely 
related In structure, de Rujula and Glashow (1980) relate the two decays by 


T « (-sr)^ t(m ♦ ey) 
m 


i"'y»3 t(m evy) 
'nT’* B(u ♦ ey) 


(37) 


which. In our discussion, corresponds to a common 6 and G'. The lower limit 
on the M ♦ ey branching ratio of 1.9 x 10”^® (Bowman, et al . 1979) and the 
p ♦ evv lifetime of 2 x 10“^ s combine to yield 


T 


1028 s 


LmMeVir 


(38) 


For chiral theories, G' (rather than G) might be the same which would 
give T » (m /m' )^ x(u ♦ ey) > 10^^/[m* (eV)]S. However, neutrinos have no 

p - 

charge, far smaller mass, perhaps Majorana character, and thus could have 
different constituents and selection rules. (Even In the standard models we 
have problems: The chiral heavy lepton model gives G'(v)/G'(p) « 0(10®).) We 


shall thus Ignore any direct p ♦ ey formula. 

Three composite model categories can be defined, showing a range of 
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lifetimes. (Thus, we cannot make a firm prediction. However, our main point 

is that composite models can give lifetimes which are much shorter than those 

predicted by the standard GWS model). The first category corresponds to 

second-order (or higher) radiative transitions where the amplitude 

< (em/A )^ and t > 10^®/[m‘ (eV)]*’s from (31) and (36). This includes models 

where the heavier quark and lepton generations are viewed as radial 

excitations (e.g., Ansel 'm (1980), note that 0(q”) is effectively 0(m''). The 

second category corresponds to 0(em/AA') amplitudes. If K' were as low as a 
30 5 

GeV, T i 10 7[m' (eV)] s. The manner in which the chiral symmetries are 

broken (and hence the size of the number c = m'/A') is very uncertain. The 

third category corresponds to models with no chiral suppression so naively 
12 3 

T > 10 /[m‘ (eV)] s. As an example, we may estimate the transition magnetic 
moment by adapting the anomalous magnetic ncment calculations of Shaw, 
Silverman and Slansky (1980) for bound states of a fermion and a boson with 
masses m^ << m^j. We find. 


^ 16 it' 


m. 

( — 2 ") (1 + tn r) 


(39) 


2 2 

where r is defined to be m^/mjj << 1. Here q and q' are the couplings between 

the neutrino composites and the two-particle states, and we have chosen 

mf « mj). If m^, * 10^ TeV, m^ = 10^ TeV and g = g' » 1 as some sort of 

22 3 

hyperstong interaction, then we get t * 10 /[m'(eV)] s. 

The smallest lifetime in those models whose neutrinos are fundamental 
fields corresponds to estimates like (28), yet still appears to be too large 
to account for any cosmic UV background flux observed. We propose here that 
significantly smaller lifetimes can be found in the case where the neutrino is 
not elementary, and that cosmic UV observations may give the first evidence 
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for composite structure of leptons. 


V. ULTRAVIOLET BACKGROUND DATA 


The observational situation regarding the cosmic ultraviolet background 
fluxes, particularly at high galactic latitudes, Is still In a relatively 
primitive state owing to fundamental observational difficulties. These 
observations have been reviewed quite recently (Paresce and Jacobsen, 1980; 
Henry 1981) and the reviews point out, among other things, conflicts In both 
observations and Interpretation. Nevertheless, the contributions from various 
sources of background contamination can be estimated and general cosmic flux 
levels can be established. Although It was originally suggested ty 
de Rujula and Glashow (1980) that the UV flux from decay of neutrinos In the 
galactic halo would have a peak Intensity In the direction of the galactic 
center, fluxes from stars and a large dust opacity make searches In this 
direction Impractical. Rather, one should look In the direction of the 
galactic poles where these effects are minimized. Indeed, significant 
portions of the sky near the galactic poles may be almost totally free of dust 
(Shane and Wirtanen 1954; Helles and Jenkins 1976) 

The UV observations (Maucherat-Joubert, Gruvelller and Deharveng 1978; 
Anderson et al . 1979) may be summarized as follows: With all nianbers In units 
of photons cm“2 s"^ sr"^ the diffuse high-latitude far UV spectrum 
appears to be flat between - 1300 X and - 1525 X with an Intensity of 260 ± 

40. (Allowance for up to 0.2 mag of extinction by high latitude dust could 
bring this number up by as much as 20 percent, but this is still within the 
error of the measurements.) In the range between 1680 ^ and 1800 the mean 
flux level Increases to ~ 600. The big question here is how much of the flux 
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could be from such things as scattered starlight, airglow, and the Integrated 
flux of distant galaxies. It has been argued that backscattering of starlight 
is negligible (Henry 1981). The ~ 1700 t feature is not consistent with 
calculations of the spectrum from distant galaxies but may be due to airglow 
(another point of contention). In the next section, we will use the "flat" 
flux level to derive a lower limit on the neutrino lifetime, and we will also 
discuss the possibility that the ~ 1700 X feature may be from neutrino decay 
(Stecker 1980) and the implications of this hypothesis. 


VI. ASTROPHYSICAL LOWER LIMITS ON t AND OTHER ASTROPHYSICAL IMPLICATIONS: 

By making use of equation (16), the measurements of discussed in 

section V can be used to place lower limits on t. The most stringent limits 

-5/2 

are obtained for the case n*l(Ij^«x ') and using the data at the 
shortest wavelengths. For this purpose, we take 

I ~ 200 cm'^ s”^ sr“^ A"^ (40) 

1250 A 

(Anderson et al . 1979). 

Most previous workers when using measurements or limits of the background 

radiation at various discrete wavelengths to obtain t(Eq) or T(m^) have erred 

in connecting these discrete points to generate a smooth function t(E^). This 

method can be quite misleading, as it fails to account for the fact that local 

0 

neutrino decay emission would occur in very narrow lines (AX ~ 1 A) at 
specific wavelengths not covered by the data set used (see Figure 3a). There 
is, however, a way to obtain a correct continuous function t(E^) by utilizing 
the fact that cosmological neutrinos produce a redshifted continuum spectrum 
given by equation (15). Figure (3b) shows the characteristic triangular 


20 


shaped spectrum obtained on a log - log x plot obtained from equation (15) 
If neutrino decay at an observation wavelength corresponding to point 0 Is 
responsible for the flux at 0 (solid triangle). However redshifted radiation 
from the decay of higher mass neutrinos can also account for the flux at 0 
(dashed triangle). The triangles are Inverted on a log flog graph (see 
Figure 3b). Adding together the limits thus obtained from flux measurements 
at several wavelengths gives a typical zig-zag limit function for T(m^) as 
indicated In Figure 3c. For this purpose, background fluxes and limits were 
obtained using the Infrared and optical data of Matllla (1976) and Dube, 

Wickes and Wilkinson (1979) and the UV as compiled and reviewed by Henry 
(1981). The resulting limit function from observational data over the whole 
frequency range of Interest (Infrared-optical-ultravlolet) Is shown in Figure 
4. The limits obtained from actual photon flux measurements correspond to the 
line labled SBp. For data coo^llatlons where the fluxes are given In units of 
F(e»*g cra‘2 s"^ sr’^ Hz"l), the Individual sections of SBp are given by the 
formula 


00 

where Is the energy corresponding to the frequency of the observation Vq|, 
and rv|. Is the Heavyside function: n^(x) » 1 for x > 0 and n^lx) « 0 for x < 0. 

0 

In the case where > 13.6 eV (the tyman limit X^ < 912 A) the decay 

photons are not generally directly observable (however, see footnote 41 
later), but the Indirect Ionizing properties of the photons can be used to 
place limits on the decqy time. This can be done by requiring that 
photoionization of high velocity clouds of neutral hydrogen (HI) near our 
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ga1a)^ not exceed observational limits (Melott and Sclama 1981). 

Utilizing the condition that the Ionization rate from v-decay photons not 
exceed the recombination rate, Melott and Sclama (1981) obtain the lower limit 


, . (« X ( T ,3/4 jm , h-l[l.(^,3/2] «) (42) 

91 2A 


where n^jj Is the density of Ionized hydrogen In cm"^, T Is temperature, d Is 
the distance of the cloud In kpc, ^ Is the angular extent of the cloud on the 
sky, n' Is In eV and N Is the number of Ionization per photon. Equation (42) 
gives a conservative lower limit on t of ~ 102**s If the clouds are at a 
distance of ~ 1 kpc. 

Another method of computing t from Ionization arguments Is to note that 

14 

the lifetime of the clouds > 10 s. In order for the clouds to exist In 
their neutral state, the Ionization rate must therefore be low enough such 
that the photons cannot eat through the cloud In a time Therefore, the 

flux from neutrino decay must satisfy 


.V 

Y 


^ ’’^cl ~ '^11 ^ 


(43) 


from which a rough limit Is obtained on the neutrino lifetime 
T z 4 X lO^^s 


(44) 


In agreement with that obtained from equation (42). 

The limits obtained from equations (42) and (43) are also shown In Figure 4. 

These limits can be compared with the limits given by equation (15). Equations 

0 

(42) and (43) only refer to the wavelength region x < 912 A which represents 
m' > 27.2 eV. The decay of lighter mass neutrinos, of course, will not produce 
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Ionizing radiation. It should be noted that if the high velocity clouds 
originate in the galactic plane, they could be continually in the process of 
"evaporating" by ionization once they leave the protection of the galactic 
disk. They can therefore start out with higher values of nt then observed. Also 
the corona of ionized plasma which would form around the neutral core of the 
cloud could significantly slow the ionization rate (Felten and Bergeron 1969). 
Both these considerations could make the limits obtained from equations (42) and 
(43) somewhat too restrictive, but we assume here thay they are "reasonable" to 
within an order of magnitude^. 


2. Recently, preprints by Sciama and Sciama and Melott have come to our 
attention which conjecture that decay of ~ 100 eV neutrinos could provide a 
source for ionizing the intergalactic medium and for ionization in the actic 
halo. Cruddace et al. (Astrophys. J., 187 , 497 (1974); have shown that some 
radiation at such wavelengths ('- 250 8 ) may be directly observable in very 
restricted regions of the sky where the hydrogen column densities are known to be 
abnormally low owing to the opacity of hydrogen dropping off steeply with energy 
for photon energies above the lyman limit. Another recent discussion of 
ionization of galactic and intergalactic gas by photons from this decay of 
30eV-150eV neutrinos is given by Raphael i and Szalay (Preprint HSF-ITP-81-52). 
They find t z 10^^ s in agreement with the results shown in Figure 4. They also 
show that in the case t « 10^^ s (higher ionizing fluxes) would have serious 
consequences for the evolution of the universe. 


Recently, Shipman and Cowsik (1981) and Henry and Feldman (1981) have 
employed another method to obtain limits on t from UV observations. These 
authors have used data on the UV flux spectrum in the directions of the Virgo and 
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Coma clusters. Assunming that this radiation may be partially due to neutrino 
decey and assunming that the (dynamical mass of these clusters Is dominated 
neutrinos of a given mass equal to twice the energy of the UV photon (see 
equation (9)), equations (11) and (13) can be used to determine a lower limit on 
T as a function of UV photon energy. The limits obtained by these authors are 
also shown In Figure 4. 

Having simim&rlzed the limits on t In Figure 4, we now discuss the 
Interesting conjecture that the ~ 1700 t feature (see section V) could be due to 
neutrino decey (Stecker 1980). This feature could then be hypothesized to be 
from a decay line somewhere In the band pass region of the photometers of 
Maucherat-Joubert, et al . (1978) and Anderson et a1 . (1979a, b) I.e., In the 
wavelength range 1680 X-1800 i! corresponding to an energy range 6. 9-7.4 eV and a 
neutrino mass m' In the range 13.8-14.8 eV. Of course, such neutrinos would have 
all of the desirable cosmological p'opertles discussed In Section II by 
satisfying the condition (6). The line would have an expected width ~ 2 ^ and 
for neutrinos In a large galactic halo would require a neutrino lifetime « 6 x 
10^** s with the point s shown on Figure 4. (Although It Is not too clear on the 
figure, this wavelength range was not covered In the observation of I'enry and 
Feldman (1951)). This lifetime Is within the limits obtained from our 
astrophysical arguments; however. It Is much shorter than that given by the 
“standard” calculations (see Figure 4). But within the framework of the new 
substructure models for leptons and quarks (see Section IV) such decay rates are 
possible, (although not required). 

Thus, If the - 1700 t feature or some similar feature, shown by future 
observations to be narrow, could be shown to be from neutrino decqy. It would be 
a test which would determine neutrino mass from equation (8) or (9) and may be 
the best way to prove that substructure for leptons and quarks or other new 
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physics exists. We therefore urge that Improved high galactic latitude searches 
be made «f1th a fleld-of-vleM small enough to exclude hot stars and dust patches 
and with good spectral resolution^. We alto suggest that such searches should 
begin with the 1680 X-1800 i region^. 


3. Conversely. Henry (1980) has suggested going to a very large (40^ diameter) 

fleld-of-vlew and subtracting out a contribution from stars of ~ 500 cm'^s’Ur’ 
^4’^. This would minimize the error from Individual faint stars. In looking for 
emission from galaxy clusters such as Virgo and Coma, one should go to a 

f1eld«of-v1ew comparable with the size of these objects. I.e.. of the order of a 

degree. 


4. Fritzsch ( Proc. Oxford Inti. Symp. on Progress In Cosmology, 1981) has 
reported that there Is now experimental evidence obtained by the CDHS accelerator 
group that the weak Intermediate vector boson has a mass greater than 100 GeV, In 
accord with predictions based on composite models tut In disagreement with the 
standard QilS p*^ed1ct1on of 83 GeV. If this result should hold up. It would 
provide a strong motivation for considerations of the astrophysical Implications 
of composite models such as those given here. 
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FIGURE CAPTIONS 

Figure 1. Nass/Luninosity ratio in solar units as a function of cosmic 

scale size. For ext»*aga1actic objects the dependences on are 
as shown on the scales. 

Figure 2. Feynman diagrams for radiative neutrino decey for GWS models with 
neutrino mixing. 

Figure 3. Improper and proper methods for obtaining t(Eq) and T(m^). 

(a) Given a discontinuous set of data points 0,0',0",... 
for I^ various X, one cannot smoothly interpolate to get 
x(mj (see text), (b) Cosmological continuum spectrum for 

redshifted emission generated by higher mass ( ) and 

minimal mass ( ) neutrinos to account for observation 0 

and resulting T(m^) limits \see text), (c) Limits obtained 
from a set cf observations 0,0',0",... using the construction 
shown in (b). 

Figure 4. Theoretical model predictions for T(m^) and astrophysical lower 

limits on hgtlEQ). (It is assumed that m^ » 2 Cq, see equation (9). 
The limits marked SBj: (Stecker-Brown, this work) were obtained 
directly from cosmic photon fluxes. The limits MSj (Melott and 
Sciama 1981) and SBj (this work) are from ionizing flux limits 
(see text). The point S is obtained from the - 1700 X feature 
(Stecker 1980). The limits marked SCC and SCV were obtained by 
Shipman and Cowsik (1981) from observations of the Coma cluster 
and the Virgo cluster. Limits obtained from other observations 
of Coma and Virgo by Henry and Feldman (1981) are labled HC and HV, 
respectively. 
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Fig. 1 









29 


R. W. Brown, Physics Department, Case Western Reserve University, 
Cleveland, OH 44106. 

F. W. Stacker, Code 665, Laboratory for High Energy Astrophysics, 
NASA/Goddard Space Flight Center, Greenbelt, MO 20771 


